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Abstract 



Stationary, axisymmetric, vacuum, solutions of Einstein's equa- 
tions are obtained as critical points of the total mass among all ax- 
isymmetric and (t, 4>) symmetric initial data with fixed angular mo- 
mentum. In this variational principle the mass is written as a positive 
definite integral over a spacelike hypersurface. It is also proved that 
if absolute minimum exists then it is equal to the absolute minimum 
of the mass among all maximal, axisymmetric, vacuum, initial data 
with fixed angular momentum. Arguments are given to support the 
conjecture that this minimum exists and is the extreme Kerr initial 
data. 

1 Introduction 

In an axisymmetric, vacuum, gravitational collapse the total angular mo- 
mentum is a conserved quantity. Therefore, if we assume, according to the 
standard picture of the gravitational collapse, that the final state will be a 
Kerr black hole the following inequality should hold for every axisymmetric, 
vacuum, asymptotically flat, complete, initial data set 





where m is the mass of the data and J the angular momentum in the asymp- 
totic region. Moreover, the equality in (0) should imply that the data is an 



slice of the extreme Kerr black hole. A counter example to {TJ will provide 
a regular vacuum data that do not settle down to a Kerr black hole. For a 
more detailed discussion of the motivations and relevance of ([TJ) and related 
inequalities see [TO], [12] an d |H]- 

Inequality ((TJ) is a property of the spacetime and not only of the data, 
since both quantities J and m are independent of the slicing. It is in fact 
a property of axisymmetric, vacuum, black holes spacetimes, because a non 
zero J (in vacuum) implies a non trivial topology on the data and this is 
expected to signal the presence of a black hole. Note, however, that the 
mass in Q is a global quantity but the angular momentum is a quasilocal 
quantity because we have assumed axial symmetry. Without axial symmetry 
we still have J defined as a global quantity at spacelike infinity, but (UJ is 
not longer true in this case. A more subtle question is whether (Q) is true 
where both m and J are quasilocal quantities, that is, whether ((TJ) is in fact a 
quasilocal property of the black hole. In general there is no unique definition 
of quasilocal mass (see the recent review on the subject [17 ). However, 
a remarkable counter example was found in [JJ in which there is a clear 
quasilocal mass definition (the Komar mass) and inequality JTJ) is violated 
at the quasilocal level. Finally, let us note that (£Q) is false for black holes in 
higher dimensions (see, for example, fH] and reference therein). 

The inequality (£Q) suggests the following variational principle: 

(i) The extreme Kerr initial data is the absolute minimum of the mass 
among all axisymmetric, vacuum, asymptotically flat and complete ini- 
tial data with fixed angular momentum. 

So far, there is no proof of ((TJ). A promising strategy to prove it is to use the 
variational formulation (i). In this article we will prove the following results, 
which are a step forward in this direction. 

The first result is the following related variational principle: 

(ii) The critical points of the mass among all the axisymmetric, (t, <p) sym- 
metric, asymptotically flat data are the stationary, axisymmetric solu- 
tions. 

A spacetime is defined to be (t, <fi) symmetric if it is symmetric under a 
simultaneous change of sign of the time coordinate t and the axial angle 
0. A data is called (t, 0) symmetric if its evolution is a (t, <fi) symmetric 
spacetime. These data are also known as "momentarily stationary data" (see 
[2] for more details). The variational principle (ii) was proved by Bardeen 
who also included matter in the formulation. It was also studied by Hawking 
[TT] for black holes including boundary terms. However, in all these works 
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the mass is not written as a positive definitive integral (see the discussion 
of section VIII in [2]). Therefore, it is not possible to relate (ii) with (i) in 
these formulations. In this article we will prove (ii) using the mass formula 
discovered by Brill [3], which is a positive definitive integral over the slice. 
Using this formulation of (ii) we will be able to prove the following: 

(i ! ) If the absolute minimum of the mass among all axisymmetric, (t, 0) 
symmetric, vacuum, asymptotically flat and complete initial data with 
fixed angular momentum exists, then it is equal to the absolute mini- 
mum of the mass among all maximal, axisymmetric, vacuum, asymp- 
totically flat and complete initial data with fixed angular momentum. 
Moreover, the absolute minimum is stationary. 

That is, we have essentially reduced the variational problem (i) to the (t, 0) 
symmetric case. Note that we have included in (i') the condition that the 
data are maximal (i.e. the trace of the second fundamental form is zero). 
This is a technical assumptions which simplifies considerably the analysis, 
but the statement is expected to be valid without it. 

There exist other variational formulations of the stationary, axisymmetric, 
equations, see [14, [16J. Particularly interesting in the present context is 
the variational formulation given by Carter [7] which is based in the Ernst 
formulation [§]. There exist a remarkable connection between (ii) in the 
present formulation and Carter's variational principle, we will prove that the 
Lagrangians differ only by a (singular) boundary term. 

2 Axially symmetric initial data and Brill proof 
of the positive mass theorem 

In this section we review Brill's positive mass theorem for axisymmetric data 
jl]. The original proof was for time-symmetric data in K 3 , here we slightly 
extend it to include maximal data and non-trivial topologies. 

An initial data set for Einstein's vacuum equations consists in a 3-manifold 
S, a Riemannian metric h a i,, and a symmetric tensor field K ab such that the 
vacuum constraint equations 



are satisfied on S; where D a and R are the Levi-Civita connection and the 
Ricci scalar associated with h ab , K = h ab K ab , and the indexes are moved 
with the metric h a b and its inverse h ab . 



D b K ab - D a K = 0, 

R + k 2 - k ab k ab = o, 



(2) 
(3) 



3 



We will assume that the initial data are axially symmetric, that is, there 
exist an axial Killing vector rj a such that 

£r)h a b = 0, £ v K ab = 0, (4) 

where £ denotes the Lie derivative. The Cauchy development of such initial 
data will be an axially symmetric spacetime. 

The Killing vector if 1 is assumed to be orthogonal with respect to h a b 
to a family of 2-surfaces in S. Under these conditions, the metric h a b can 
be characterized by two functions, one is essentially the norm of the Killing 
vector and the other is a conformal factor on the 2-surfaces. We make explicit 
the choice of the free functions as follows. Let (p, z, 0) be local coordinates 
in S such that the metric has following form 

Kb = ?p A h ab , (5) 

where the conformal metric h ab is given 

h = e - 2q (dp 2 + dz 2 ) + p 2 V, (6) 

and q, ip are functions which depend only on z and p with tp > 0. The 
vector 7] a = (d/dip) a is a Killing vector of both metrics h ab and h ab . The 
norm of rf with respect to the physical metric will be denoted by X, (i.e. 
X = i] a r] b h a b = ip 4 p 2 ), the norm of r) a with respect to the conformal metric 
is given by p 2 = rj a ifh ab . 

We define the following quantity 

J(E) = £n ab r] a n b ds- h , (7) 

where 7r a f, = K ab — h a bK, S is any closed 2-surface, h a is the unit normal 
vector to S with respect to h a b and ds^ is the area element of S with respect 
to h ab . Equation (J2J and the Killing equation imply that the vector 7t a bV a 
is divergence free. If S is the boundary of some compact domain Q C S, 
by the Gauss theorem, we have J(E) = 0. For example, if S = M 3 then 
J(S) = for all S. In an asymptotically flat data, J(Soo) gives the total 
angular momentum, where Soo is any closed surface in the asymptotic region. 
Then, the angular momentum will be zero unless is not the boundary of 
some compact domain contained in S. 

In order to have non zero angular momentum we will allow S to have many 
asymptotic ends 1 . Let ik a finite number of points in R 3 . The manifold S is 

1 There is an interesting alternative (not included here) discussed in ^D] to allow non 
zero angular momentum: the interior of the manifold is assumed to be compact and non 
simply connected with a pseudo axial Killing vector. 
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assumed to be R 3 \ J2k The points i>k will represent the extra asymptotic 
ends, at those points we will impose singular boundary conditions for ip. The 
be consistent with the axial symmetry assumption the points if. should be 
located on the axis p = 0. 

In addition to axial symmetry we will assume that the data are maximal 

K = 0. (8) 

By equation (J3J this implies that R is positive, this will be essential in order 
to extend Brill's proof to non-time symmetric data. 

Define the conformal second fundamental form by K ab = ip 10 K ab . Using 
(JEJ) and flU) we obtain 

£ v K ab = 0, K = 0. (9) 

The constraint equations (0)-© can be written as equations for K ab and 
ip using the well known conformal method (see, for example, [2] and reference 
therein) 

D a K ab = 0, (10) 
D a D a il> - \rP = - X -K ah K ah ^-\ (11) 

o o 

where D a and R are the Levi-Civita connection and the Ricci scalar asso- 
ciated with the conformal metric h^. In these equations, the indexes are 
moved with the conformal metric and its inverse h ab . 

The function q is assumed to be smooth with respect to the coordinates 
(p, z). At the axis we impose the regularity condition 

q(p = 0,z) = 0. (12) 

Note that condition (fT2"|) includes the points ik- These points are assumed to 
be regular points of the conformal metric h a b, that is, h a b is well defined in 
R 3 . 

We assume the following fall-off condition at infinity 

q = o(r _1 ), q. x = o(r~ 2 ), (13) 

where r = a/ p 2 + z 2 and a comma denotes partial derivatives. This fall off 

conditions imply that the total mass of the conformal metric is zero. 

At infinity, the conformal factor ip and the conformal second fundamental 
form satisfy 

tf;-l = 0(r~ 1 ), ijj ir = 0(r~ 2 ), (14) 

and 

K ab = 0{r- 2 ). (15) 
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Under these assumptions the total mass of the physical metric is given by 

m = — — lim (p n a D a tp dsh, (16) 



2tt 



r—too 



where S r are the 2-spheres r = constant, n a is the unit normal, with respect 
to h ab , pointed outwards and dsh is the area element of E with respect to 
h a b- 

The Ricci scalar R of the conformal metric (jOJ) is given by 

R = 2e 2q {q, pp + q. zz ). (17) 
We have the important equation 

Rdp h = 0, (18) 



where dph is the volume element of the metric h a b To prove this, note that 
dfjLh = pe~ 2g dpdzdcj), then 

p roo /*oo 

I Rdp h = 4tt / dpi dz(q :PP + q tZZ )p (19) 
Jr 3 Jo J-oo 

/*oo /*oo 

= 4vr / dp dz ((pq, P -q), P + (pq,z),z), (20) 

JO J-oo 

we use the divergence theorem in two dimension to transform this volume 
integral in a boundary integral over the axis p = and infinity. The bound- 
ary integral at the axis vanishes since q satisfies ()12j) and at infinity it also 
vanishes because of (JT3J). 

Since lim,.-,.^ ^ = 1, we have an equivalent expression for the mass 



-1 r n^Dal 
m = — hm (h - ds h . (21) 

27T r^oo J lp 



We use the identity 



a v ip )~ v r ' 1 ] 

the constraint equation equation ()18j) and the mass formula (|21J) to 

obtain the final expression 

m =27L(— + —^-)^ (23) 
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which is definite positive. To obtain ()23|) from (}2*2~j) we have assumed that 
the boundary integral around the singular points vanishes, that is 

lim <b ds h = 0, 24 

where is the distance to the point ik- This condition (which is, of course, 
trivially satisfied when the topology of the physical data is M 3 ) allows for a 
singular behavior of ib at ik which in particular include the case where ik are 
asymptotically flat ends. Near an asymptotically flat end if. the conformal 
factor satisfies ib = O^ 1 ), ip^ r = 0(r^ 2 ) which imply pijl . To illustrate 
this, consider the following two examples. 

The Schwarzschild initial data in isotropic coordinates is time-symmetric 
(K a b = 0) and conformally flat (q = 0). In this case we have one point io 
located at the origin and the conformal factor is given by 

* = ! + (25) 

where is the Schwarzschild mass. We have 

f°° (it ) 2 

m = 2 ±-^-r 2 dr, (26) 
Jo V 

= m . (27) 

Note that the integral is taken over the two asymptotic regions. 

The second example is the Brill-Lindquist[5 initial data. In this case the 
data is also time-symmetric and conformally flat, but here we have n ends ik 
and the conformal factor is given by 

n 

k 

where rrik are arbitrary positive constants. The conformal factor (j28|) satisfies 
()24j) and we have that 

1 f ( Dg^D^ \ A 

m =toh{—^) d » h = 2-< mk - (29) 

In the non time-symmetric case, we have assumed that the integral of 
K ab K a i)ib~ 8 over M. 3 is bounded. At infinity, the integral converges because 
the assumptions (|TK|) and ((HJ). At the points ik the conformal second 
fundamental form will, in general, be singular. However the integral will 
be bounded because the singular behavior of K ab will be canceled out by 
the singular behavior of ib. For example, in the asymptotically flat case, 
K a b = 0(r^ 4 ) near ik and then we have that K ab K a b^~ 8 is bounded. In 
appendix El we prove that Kerr initial data satisfy these conditions. 
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3 The variational principle 



In the integral (f2*3*|) the mass depends on the metric variables ip, q (the 
function q appears in the volume element and in the indexes contractions) 
and on the conformal second fundamental form K ab . These functions are not 
independent, they have to satisfy the constraint equations (fTUj) and (fTTj) . In 
order to formulate the variational principle we want to express the mass in 
terms of functions that can be freely varied. We analyze first the conformal 
second fundamental form K ab and the constraint 
Consider the following vector field S a 

S a = K abV b - p-\ a K h 4rt. (30) 

Using equations (jHJ), (fTIl|) and the Killing equation for rf it follows that S a 
satisfies 

£^ = 0, S a Va = 0, D a S a = 0. (31) 
From (J7|) we deduce an equivalent expression for the total angular momentum 

J = / S a n a ds h} (32) 

where we have used that the second term in the right-hand side of ()30|) does 
not contribute to the angular momentum because we can always chose a 
closed surface at infinity such that n a r] a = 0. 

The conformal metric h ab can be decomposed into 

hob = q a b + P~ 2 VaVb, (33) 

where 

q ab = e ~ 2q (dp 2 + dz 2 ), (34) 

is the intrinsic metric of the planes orthogonal to r] a . Using this decomposi- 
tion and the definition of S a we obtain the following expression for the square 
of the conformal second fundamental form 

K ab K ab = K ab K cd q ac q bd + p-\K ab r, a r, b ) 2 + 2p~ 2 S a S a . (35) 

The two first terms in the right hand side of this equation are positive, then 
we have 

K ab K ah > 2p- 2 S a S a . (36) 

Equations ()32|) and (|36j) are important because they show that S a contains 
the angular momentum of K ab and its square is a lower bound for the square 
of K ab . 
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We define the tensor 

R ab = -S {a r] b \ (37) 
V 



we have 



K ab K ah = (3* 
P 



It is interesting to note (but we will not make use of it) that this tensor is 
trace free and divergence free. To prove this we use the Killing equation 
D(aVb) = 0, the fact that rf is hypersurface orthogonal, (i.e.; it satisfies 
D a i]b = — V[aDb] In rj) and equations (I3*B . 

A data will be (t, 0) symmetric if and only if the following conditions hold 
(see 0) 

K ab q ac q bd = 0, iWV = 0. (39) 

This is equivalent to K ab = K ab . 

The vector S a can be expressed in terms of a free potential. Define the 
rescaled vector s a by 

then 

£ v s a = 0, s\ = 0, d a s a = 0, (41) 
where d a is the connexion with respect to the flat metric 

5 = dp 2 + dz 2 + p 2 dif 2 , (42) 

and in equation (}4~Tj) the indexes are moved with this metric and its inverse. 
The same will apply to all the equations from now on: all of them will be 
given in term of the flat metric 5 a b and its connexion d a . 

An arbitrary vector s a , which satisfies equations (j41j) . can be written in 
term of a potential Y in the following form 

s a = ^abc^g^ (43) 

2p z 



where e a b c is the volume element of the flat metric ()42|) and £ V Y = 0. The 
motivation of the normalization factor 1/2 in will be clear in the next 
section. We have the relation 



ra^H^aw (44) 



The angular momentum (}3*Tj) is given in terms of the potential Y by 

J=l{Y(p = Q,-z)-Y(p = 0,z)), (45) 
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where z is taken to be larger than the location of any point i^. 

Motivated by Brill's formula (j2HJ), we define the mass functional as follows 

M{v, Y) = J- [ (I6d a vd a v + p- 4 e- Sv d a Yd a Y) dp , (46) 
62-k J R3 

where v = ln^ and dpio is the flat volume element. Note that in the integral 
(|4l)|) the metric function q does not appear. 

From equation (}23|) and (jH)) we see that for every axisymmetric and 
(t,<p) symmetric data we have m = Ai(v,Y). From (fHfij) we see that for 
every axisymmetric, maximal data, we have 

m>M{v,Y). (47) 

We emphasize that the functions (v , Y) can be computed for an arbitrary ax- 
isymmetric data (in the construction of the potential Y we have not used the 
maximal condition) and then the functional Ai(v,Y) can be also calculated 
for arbitrary data (provided, of course, the integral is well defined). However, 
only for maximal data we can use the Brill formula (f2*3|) to conclude (|47|) and 
only for (t, 0) symmetric data we have that M.(v , Y) is in fact the mass. 

For the present calculations is more convenient to write the functional Ai 
in the form ([46)1 . where the axial symmetry is not explicit. For completeness, 
we also write it in a manifest axisymmetric form 

-| /»oo poo 

M(v,Y) = -J dp J dz(lQp(v 2 z + vl)+p- 3 e-^{Y" z + Y 2 p )). (48) 

Let us define A as the set of all functions (v, Y) such that the integral (f4T)J) 
is bounded. Although A4(v, Y) is well defined in A, not for every function in 
A we will have that M.(v, Y) is equal to the mass of some (t, <p) symmetric 
initial data. This is a subtle and important point, let us discuss it in de- 
tail. We have seen that all axisymmetric and (t, 0) symmetric data can be 
generated by three functions (v , q, Y). They are coupled by the Hamiltonian 
constraint (JSJ). In coordinates, this equation is given by 

Aip . . d a Yd a Y . . 

where A is the flat Laplacian with respect to (|42jl . For given (v, Y) (remember 
that v = \nip) this is a linear, two dimensional, Poisson equation for q. The 
delicate point are the boundary conditions. In order to obtain Brill's formula 
we have required that q satisfies (JT2J) and (fT3|). But we cannot impose this 
two equations as boundary conditions for a two dimensional Poisson equation. 
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Let say that we impose ()12j) and we ask for solutions which fall off at infinity. 
This problem can be solved with an explicit Green function. However, in 
general, the fall off of the solution will be q = 0(r _1 ) which is weaker than 
(JTHJ). Only for some particular source functions {v , Y) the solution q will 
satisfy (JTHJ). Let us denote by A\ the subset of A of those functions {y , Y) 
such that the solution q of equation (pf§j) satisfies (fT2j) and (|T3*j) . Only for 
functions in A\ the functional A4(v,Y) can be written as a the boundary 
integral (fT6j) and hence gives the mass of some initial data. A function v of 
compact support (such that ip = I near infinity) is an example of a function 
which is in A but not in Ai (we can take Y — 0), since in this case clearly 
A4(v,Y) is strictly positive and the boundary integral (fT^j) is zero. 

We want to make variations of Ai(v,Y). At first sight, it appears that 
the appropriate set for admissible functions is A± and not A. However, it 
seems to be difficult to characterize A\. It is known how to characterize the 
set of those q such that (jUIJ) has a solution if) (for an, essentially, arbitrary 
Y) which satisfies dEJ), in this case a non-linear equation must be solved (see 
jE] and PH]). However, this set is not very useful in the present context since 
for the Brill formula is natural to use (v , Y) as independent functions and 
not (q, Y). Instead, what we will do is to take A as the set of admissible 
functions. Remarkably, it will turn out that the critical equations in this 
bigger set are only the stationary, axially symmetric equations. 

Let a and y be compact supported functions in M 3 with support in S and 
such that the support of y does not contain the axis. By equation P5^) we 
see that this condition implies that the perturbation Y + y does not change 
the angular momentum of Y. Define 

i(e) = M(v + ea,Y + ey). (50) 

The first variation of Ai(v , Y) is given by 

i'(0) — — ~~~ — f (I6d a vd a a - Aap- 4 e - 8v d a Yd a Y + p- 4 e - 8v d a Yd a y) dfi , 

167T Jr3 

(51) 

where a prime denotes derivative with respect to e. Integrating by parts, we 
obtain that the condition 

i'(0) = 0, (52) 
for all a and y is equivalent to the following Euler-Lagrange equations 

AAv + p- 4 e~ 8v d a Yd a Y = 0, (53) 
d a (p- 4 e~ 8v d a Y) = 0. (54) 
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The second variation is given by 

i"(0) = ~~~ f {16d a ad a a+ 
16vr J r3 

(32a 2 d a Yd a Y - lQad a Yd a y + d a yd a y) p~ A e - Sv } dfx - (55) 

There is an equivalent way of deducing equations (JH3J) - (JH4*|) . Instead of 
taking Y as variable we take the vector s a , which should satisfy the con- 
straints (j41j) . The mass functional is given by 



M(v, s) = ±- [ {Ad a vd a v + p- 2 e- 8v s a s a ) dp . 



(56) 



Let 7° be a compact supported vector in S such that the support of 7° does 
not contain the axis. We assume that 7° satisfies the constraint 

dal a = 0. (57) 
We define i in analogous way as in (|50j) . The first variation is given by 



i'(0) = ±- ! {Ad a vd a a - Aap~ 2 e' Sv s a s a + p- 2 e^ v s a la ) dp , 
integrating by parts we get 

= 7- f (-4a(A« + p- 2 e- 8 ^ a s a ) + /rV 8 ^) <W 
From this we deduce the Euler-Lagrange equations 



(5f 



(59) 



Av + p- 2 e-* v s a s a = 0, (60) 

p- 2 e-*>s a = ±d a Q, (61) 

for some function Q. Equation ((HOj) follows because we can make arbitrary 
variations in a. On the other hand, variations in 7° should satisfy the con- 
straint (|57j) . Writing 7° as the curl of an arbitrary vector and integrating by 
parts we get 

d [a H b] = 0, (62) 

where 

H a = p- 2 s a e- 8v . (63) 

Equation (Jo^j) is equivalent to (JBTj) . Using the constraint <9 a s a = 0, we deduce 
the following equations which does not involve s a 

4Av + p 2 e Sv d a Qd a Q = 0, (64) 
d a (p 2 e Sv d a n) = 0. (65) 
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Equations (jo^|) - (|63j) are equivalent to equations (px?j) - (p>^) . the relation be- 
tween Q and Y is given by 

d a Q = p- 4 e- Sv e abcV b d c Y. (66) 

In the next section we will prove that these equations are precisely the sta- 
tionary, axisymmetric, vacuum equations. This will provide also an interpre- 
tation for the potential Y and the velocity Q in the stationary case. Note 
that Y is defined for arbitrary data, in contrast Q is only defined for solutions 
of the critical equations, that is, for stationary axisymmetric data. 
If we take Y = 0, then these equations reduce to 

Av = 0, (67) 

which is Weyl equation for axisymmetric, static, spacetimes. This is of course 
consistent with the result that we are going to prove in next section. However, 
it is important to note that the Schwarzschild data in the form (|25|) does not 
satisfy (|6T|). Schwarzschild satisfies (|57jl in Weyl coordinates where v and 
the metric function q are given by 

1, / f, + f _ — 2m\ 1 f (r+ + f) 2 — Am 2 
v = —— In — = , q = - In — — 

4 \r + + r_ + 2m J 2 \ Ar + r_ 

with r± = p 2 + (z ± m) 2 . The relation with the isotropic coordinates (r, 9) 
used in (|23j) is 

(77l 2 \ ( TTL 2 \ 

l -\T 2 )> ~ Z = Z { 1 + IT 2 )> m 

where z = rcosO and p = rsin^. Since X is an scalar independent of coor- 
dinates we have X = p 2 tp A = p 2 tp 4 . The function q satisfies our assumptions 
(p"2|) and (|T3|) . however the conformal factor ip = e v does not satisfies (j2lj) . 
The conformal factor is singular on the rod p = 0, — m < z < m (which 
represent the horizon of Schwarzschild data) and not just on singular points 
ik- The integral A4(v,0) diverges. Note that 1R 3 in Weyl coordinates (p,z) 
represent the exterior of the black holes, in contrast to coordinates (p, z) 
where M 3 represent both asymptotic regions. 

4 Stationary axisymmetric fields 

The spacetime metric of a vacuum, stationary and axially symmetric space- 
time can be written, in Weyl coordinates, in the following form (see, for 
example, fTSI) 

g = -V(dt - adcf)) 2 + V- 1 [p 2 d(f) 2 + e 2 ^(dp 2 + dz 2 )] , (70) 
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where the functions V, a and 7 depend only on (p,z). The two Killing 
vectors are 

'-(£)'• 

they define the scalars 

V = -eC9,u, X = rfrfg^ W = rf^g^ (71) 

where /1, v are spacetime indexes. We have the following relations 

W = Va, p 2 = VX + W 2 . (72) 

The vacuum field equations are given by 

d a {V- l d a V + P ~ 2 V 2 ad a a) = 0, (73) 
d a {p- 2 V 2 d a a) = 0. (74) 

We want to prove that these equations are equivalent to equations (Jfi4]) - (fB5j) . 
We first compute the relation between (V, a) and (v, fl). 

Take an slice t = constant of the metric (|70j) . The intrinsic metric of this 
surface is given by 

h = V-\p 2 - a 2 V 2 )d<f) 2 + V~ x c 2 \dp 2 + dz 2 ). (75) 

To write this metric in to the form ([6]) set 

r= (P 2 - V2 ^ =X= X ( 76 ) 

v v P 2 p 2 {vx + w 2 y { } 

and 

e 2 i = —P- = tJL (77) 

{p 2 -V 2 a 2 ) VX ' { J 



From (fTHJ) we deduce 



In order to compute Q(v, a) we need to calculate the second fundamental 
form of this foliation. The lapse and the shift of the foliation t = constant 
are given by 

N = -£= = V" 2 , N a = W(d<P) a . (79) 
v X 
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and the second fundamental form is 

K ab = -^D {a N b) . (80) 

We write N a in terms of the Killing vector r] a , as in the previous section we 
define fj a = h ab rj b where h ab is given by (J7H|) . then we have 

Va = ^— ^ l {d<j>) a . (81) 

Using this expression we write N a as 

N a = Qfj a , (82) 

where Q is given by 

V 2 a W 

^" ^-^ T (83) 

The scalar can be interpreted as the angular velocity of the locally non 
rotating observers (see [2] and also [TH] p. 187). We want to prove that this 
function Q is precisely the potential Q of the previous section. In order to 
see this let us compute the vector s a 

s a = V b K ab = ^ 2 V b K ab = ^ 2 Xd a Q = l -^p 2 d a n. (84) 

Where we have used s a = s b 5 ab = S b h ab . Equation (Jgljl is identical to 
equation (foTj). 

Using the relations (J7%|) and (jHSJ), after a long but straightforward com- 
putation, we conclude that equations (J73|) and (|71jl for the functions (V, a) 
are equivalent to equations (fBO jl -lffiT j) for (v, Q). 

There is another way to prove the equivalence with the stationary equa- 
tions, using the potential Y . We replace v by X, that is we consider X, Y as 
variables. From equation ([THjl we get 

v = -InX - -In p. (85) 



Take the functional Ai defined in ()46|) but let us perform the integral on a 
bounded domain B, in terms of the variables X, Y we get 

M{X,Y)=M'(X,Y)-±- J ln(-|) Alnpd/io + y p" 1 In (£) n a «9 a pd S , 

(86) 
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where we have defined 



M'(X,Y) 



32tt 



1 



/. 



( 



daXcPX + d a Yd a Y 
X~ 2 



dp . 



(87) 



But we have 



A In p = 



8. 



.8 



for p ^ 0. Then M. and .M' differ only by a boundary term. Hence they 
give the same Euler-Lagrange equations. Note, however, that the boundary 
term is singular at the axis p = 0: if we take a cylinder p = constant near 
the axis we have, X = 0(p 2 ), ds = pdzdcj), n a d a p = 1, then the boundary 
term diverges like O(lnp) as p — > 0. 

In [7J Carter formulates a variational principle for the axisymmetric, sta- 
tionary equations. This formulation is, essentially, a modification of the [Oj 
formulation in which the norm of the axial Killing vector (and not of the 
stationary one) is taken to be the principal variable. 

Carter's Lagrangian is precisely M 1 (we use the same notation for X and 
Y, this is the reason for the normalization factor 1/2 in ([43))) . In [7J it is 
proved that the critical equations of Ai' are the stationary, axisymmetric 
equations. Therefore, the same is valid for A4. There are, however, some 
important points that we want to stress. 

If we ignore boundary terms, then equation (|86|) provides an interpre- 
tation for Carter Lagrangian. Also, it gives an interpretation of the space 
of admissible functions in which the variations are made for the following 
reason. In Carter's formulation Y is defined in terms of W and X by 



This equation can easily be obtained from (JUS!), (J7J3) and (JHBJ). That is, Y 
is defined only for stationary axisymmetric spacetimes. From the discussion 
of section El we have seen that Y can be defined for arbitrary, axisymmetric 
data, and the variation of Y and X are in fact variation among axisymmetric 
and (t, <fi) symmetric data. 

Let us consider boundary terms. The behavior of X near the axis implies 
that M 1 is singular if the domain of integration includes the axis. On the 
other hand we have seen that M. is finite. In particular, in appendix 1X1 we 
have explicitly checked that Kerr initial data in quasi-isotropic coordinates 
satisfy all our assumptions and then M. is finite and equal to the mass for 
Kerr. However, is important to note that the relevant domains of integration 
are different in Carter's formulation and in the present one. In [7J, the 
domain is the black hole exterior region, in which the inner boundary is the 
horizon. In section El we have not included any inner boundary conditions, 



e abc V b d c Y = Xd a W - Wd a X. 



(89) 
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the domain of integration is the whole manifold which can include many 
asymptotic ends. This difference is reflected in the choice of the coordinate 
system. We have discussed this with Schwarzschild data in section El The 
same apply to non-extreme Kerr initial data in Weyl coordinates: M. is 
singular in this coordinates. However, for extreme Kerr, the Weyl coordinates 
and the quasi-isotropic coordinates coincides. In this case both domains of 
integration coincides and Ai is finite whether A4' is not. 

5 Final comments 

We have analyzed the first variation of the, positive definite, mass functional 
M. (defined by (j46j) ) over axisymmetric and (t, <ft) symmetric initial data 
with fixed angular momentum. We have shown that the critical points are 
the stationary, axial symmetric equations. This proves the variational prin- 
ciple (ii). The functional is a lower bound for the mass (inequality (@Zj)) 
for all maximal, axisymmetric data. This proves (i'). In order to prove (i), 
and hence inequality we should prove that extreme Kerr is the unique 
absolute minimum of Ai over axisymmetric and (t, 0) symmetric initial data 
with fixed angular momentum. This will require the study of the second 
variation of M., given in equation 
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A Kerr initial data 

Consider the Kerr metric in Boyer-Lindquist coordinates (t,r,9,4>). The 
scalars (|TT|) are given by 



V 



A -a 2 sin 2 9 



W 



2mar sin 



2 9 



(90) 




(91) 



where 



A = r 2 + a 2 — 2mr, S = r 2 + a 2 cos 2 9, 



(92) 
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and m is the total mass and a is the angular momentum per unit mass (i.e. 
J = ma). 

The intrinsic metric hah of a hypersurface t = constant in these coordi- 
nates is given by 

h = ^dr 2 + Ed9 2 + nd(f> 2 . (93) 

The metric (193)1 has a coordinate singularity when A = 0. The solutions of 
the equation A = are given by 



+ \J m? — a 2 , r_ = m — \J m 2 — a 2 . (94) 



By the following coordinate transformation we extend the metric to a com- 
plete manifold with two asymptotic ends. Let us define the quasi-isotropic 
radius r as the positive root of the following equation 

m 2 - a 2 . 

r = r + m H . 95 

4r 

Note that when a = this reduce to the isotropic radius for the Schwarzschild 
metric. The manifold (like in the Schwarzschild case) has to isometric asymp- 
totically flat components (the region f > r + of the metric (|93j)) joined at the 
minimal surface (the horizon) r — r + . The components of h a b in the coordi- 
nates (r, 0, 0) are given by 



h = ^dr 2 + Y,dd 2 + ndcf) 2 . (96) 



The metric ()96)1 has the form (|5*|)-(|6]) with 

(9?) 

where p = rsin# and z = rcosd. Assume m > \a\. Then, from (|9Tj) we see 
that in the limit r — > we have 



y/m 2 — a 2 m . , n , 9 , 

V>=- + +O r , ^ r = O r~ 2 9£ 

r 2Vm 2 - a 2 



and at infinity 
From ()97|) we also have that 



777 

= 1+ +0(r- 2 ), g = 0(r- 2 ). (99) 
2r 



q(p = 0) = 0. (100) 
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Hence, q satisfies (JT2J), (fT3*|) and ip satisfies (JHJ) and 

The velocity fl can be calculated from equation (jH^j) using (J§Uj) and ()9"T]) 



2mar 

(r 2 + a 2 ) 2 - Aa 2 sin 2 6' 



" = - 75 ■ _^ ^—27 - (101) 



The potential K is given by 

, „ „. 2ma 3 cos 9 sin # 
K = 2ma(cos ^ — 3 cos 9) . (102) 

Zj 

Note that equation (j43j) is satisfied for z ^ 0. To see that the integral of 
d a Y d a Y p~ A ip~ 8 over IR 3 is bounded we need to check the behavior of this 
function at infinity and at the axis p = 0. At infinity we have 

d a Yd a Yp' 4 = 0{r~ 6 ), (103) 

and at the axis 

d a Yd a Yp- A ip- 8 = 0(r 2 ), (104) 

where we have used the ()98|) . Then, the integral is bounded and therefore 
we have proved that the Kerr initial data satisfies our assumptions which 
implies that M.(v, Y) = m. 

Weyl coordinates (p, z) are related to the coordinates (r, 9) by 

p = V / Asin^, z=(r — m) cos 9. (105) 

Consider now the extreme case m = \a\. In this case we have 

r = f — m, A = r 2 , (106) 

and the coordinates (r, 9) are equal to the Weyl coordinates. Equations ()100|) 
and ()99|) are still valid in this case. The fall off of the conformal factor near 
r = is however different 

this is because r = is not an asymptotically flat in this case. Nevertheless 
t/> satisfies (124)1 . The fall of behavior of Y at infinity is the same as in the 
non-extreme case. Near the axis, because of ([107J1 . we have 

d a Yd a Yp- 4 4>- 8 = 0(r~ 2 ), (108) 

and hence we conclude that d a Yd a Yp~ A ifj~ 8 is integrable over IR 3 . 
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